3.2.85 \(\int \coth ^4(x) \sqrt {a+b \text {sech}^2(x)} \, dx\) [185]
Optimal. Leaf size=84 \[ \sqrt {a} \tanh ^{-1}\left (\frac {\sqrt {a} \tanh (x)}{\sqrt {a+b-b \tanh ^2(x)}}\right )-\frac {(3 a+2 b) \coth (x) \sqrt {a+b-b \tanh ^2(x)}}{3 (a+b)}-\frac {1}{3} \coth ^3(x) \sqrt {a+b-b \tanh ^2(x)} \]
[Out]
arctanh(a^(1/2)*tanh(x)/(a+b-b*tanh(x)^2)^(1/2))*a^(1/2)-1/3*(3*a+2*b)*coth(x)*(a+b-b*tanh(x)^2)^(1/2)/(a+b)-1
/3*coth(x)^3*(a+b-b*tanh(x)^2)^(1/2)
________________________________________________________________________________________
Rubi [A]
time = 0.17, antiderivative size = 84, normalized size of antiderivative = 1.00, number of steps
used = 7, number of rules used = 7, integrand size = 17, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.412, Rules used = {4226, 2000,
486, 597, 12, 385, 212} \begin {gather*} \sqrt {a} \tanh ^{-1}\left (\frac {\sqrt {a} \tanh (x)}{\sqrt {a-b \tanh ^2(x)+b}}\right )-\frac {1}{3} \coth ^3(x) \sqrt {a-b \tanh ^2(x)+b}-\frac {(3 a+2 b) \coth (x) \sqrt {a-b \tanh ^2(x)+b}}{3 (a+b)} \end {gather*}
Antiderivative was successfully verified.
[In]
Int[Coth[x]^4*Sqrt[a + b*Sech[x]^2],x]
[Out]
Sqrt[a]*ArcTanh[(Sqrt[a]*Tanh[x])/Sqrt[a + b - b*Tanh[x]^2]] - ((3*a + 2*b)*Coth[x]*Sqrt[a + b - b*Tanh[x]^2])
/(3*(a + b)) - (Coth[x]^3*Sqrt[a + b - b*Tanh[x]^2])/3
Rule 12
Int[(a_)*(u_), x_Symbol] :> Dist[a, Int[u, x], x] /; FreeQ[a, x] && !MatchQ[u, (b_)*(v_) /; FreeQ[b, x]]
Rule 212
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
/; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])
Rule 385
Int[((a_) + (b_.)*(x_)^(n_))^(p_)/((c_) + (d_.)*(x_)^(n_)), x_Symbol] :> Subst[Int[1/(c - (b*c - a*d)*x^n), x]
, x, x/(a + b*x^n)^(1/n)] /; FreeQ[{a, b, c, d}, x] && NeQ[b*c - a*d, 0] && EqQ[n*p + 1, 0] && IntegerQ[n]
Rule 486
Int[((e_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_))^(q_), x_Symbol] :> Simp[(e*x)^(m
+ 1)*(a + b*x^n)^(p + 1)*((c + d*x^n)^q/(a*e*(m + 1))), x] - Dist[1/(a*e^n*(m + 1)), Int[(e*x)^(m + n)*(a + b*
x^n)^p*(c + d*x^n)^(q - 1)*Simp[c*b*(m + 1) + n*(b*c*(p + 1) + a*d*q) + d*(b*(m + 1) + b*n*(p + q + 1))*x^n, x
], x], x] /; FreeQ[{a, b, c, d, e, p}, x] && NeQ[b*c - a*d, 0] && IGtQ[n, 0] && LtQ[0, q, 1] && LtQ[m, -1] &&
IntBinomialQ[a, b, c, d, e, m, n, p, q, x]
Rule 597
Int[((g_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^(n_))^(p_.)*((c_) + (d_.)*(x_)^(n_))^(q_.)*((e_) + (f_.)*(x_)^(n_)),
x_Symbol] :> Simp[e*(g*x)^(m + 1)*(a + b*x^n)^(p + 1)*((c + d*x^n)^(q + 1)/(a*c*g*(m + 1))), x] + Dist[1/(a*c*
g^n*(m + 1)), Int[(g*x)^(m + n)*(a + b*x^n)^p*(c + d*x^n)^q*Simp[a*f*c*(m + 1) - e*(b*c + a*d)*(m + n + 1) - e
*n*(b*c*p + a*d*q) - b*e*d*(m + n*(p + q + 2) + 1)*x^n, x], x], x] /; FreeQ[{a, b, c, d, e, f, g, p, q}, x] &&
IGtQ[n, 0] && LtQ[m, -1]
Rule 2000
Int[(u_)^(p_.)*(v_)^(q_.)*((e_.)*(x_))^(m_.), x_Symbol] :> Int[(e*x)^m*ExpandToSum[u, x]^p*ExpandToSum[v, x]^q
, x] /; FreeQ[{e, m, p, q}, x] && BinomialQ[{u, v}, x] && EqQ[BinomialDegree[u, x] - BinomialDegree[v, x], 0]
&& !BinomialMatchQ[{u, v}, x]
Rule 4226
Int[((a_) + (b_.)*sec[(e_.) + (f_.)*(x_)]^(n_))^(p_.)*((d_.)*tan[(e_.) + (f_.)*(x_)])^(m_), x_Symbol] :> With[
{ff = FreeFactors[Tan[e + f*x], x]}, Dist[ff/f, Subst[Int[(d*ff*x)^m*((a + b*(1 + ff^2*x^2)^(n/2))^p/(1 + ff^2
*x^2)), x], x, Tan[e + f*x]/ff], x]] /; FreeQ[{a, b, d, e, f, m, p}, x] && IntegerQ[n/2] && (IntegerQ[m/2] ||
EqQ[n, 2])
Rubi steps
\begin {align*} \int \coth ^4(x) \sqrt {a+b \text {sech}^2(x)} \, dx &=\text {Subst}\left (\int \frac {\sqrt {a+b \left (1-x^2\right )}}{x^4 \left (1-x^2\right )} \, dx,x,\tanh (x)\right )\\ &=\text {Subst}\left (\int \frac {\sqrt {a+b-b x^2}}{x^4 \left (1-x^2\right )} \, dx,x,\tanh (x)\right )\\ &=-\frac {1}{3} \coth ^3(x) \sqrt {a+b-b \tanh ^2(x)}+\frac {1}{3} \text {Subst}\left (\int \frac {3 a+2 b-2 b x^2}{x^2 \left (1-x^2\right ) \sqrt {a+b-b x^2}} \, dx,x,\tanh (x)\right )\\ &=-\frac {(3 a+2 b) \coth (x) \sqrt {a+b-b \tanh ^2(x)}}{3 (a+b)}-\frac {1}{3} \coth ^3(x) \sqrt {a+b-b \tanh ^2(x)}-\frac {\text {Subst}\left (\int -\frac {3 a (a+b)}{\left (1-x^2\right ) \sqrt {a+b-b x^2}} \, dx,x,\tanh (x)\right )}{3 (a+b)}\\ &=-\frac {(3 a+2 b) \coth (x) \sqrt {a+b-b \tanh ^2(x)}}{3 (a+b)}-\frac {1}{3} \coth ^3(x) \sqrt {a+b-b \tanh ^2(x)}+a \text {Subst}\left (\int \frac {1}{\left (1-x^2\right ) \sqrt {a+b-b x^2}} \, dx,x,\tanh (x)\right )\\ &=-\frac {(3 a+2 b) \coth (x) \sqrt {a+b-b \tanh ^2(x)}}{3 (a+b)}-\frac {1}{3} \coth ^3(x) \sqrt {a+b-b \tanh ^2(x)}+a \text {Subst}\left (\int \frac {1}{1-a x^2} \, dx,x,\frac {\tanh (x)}{\sqrt {a+b-b \tanh ^2(x)}}\right )\\ &=\sqrt {a} \tanh ^{-1}\left (\frac {\sqrt {a} \tanh (x)}{\sqrt {a+b-b \tanh ^2(x)}}\right )-\frac {(3 a+2 b) \coth (x) \sqrt {a+b-b \tanh ^2(x)}}{3 (a+b)}-\frac {1}{3} \coth ^3(x) \sqrt {a+b-b \tanh ^2(x)}\\ \end {align*}
________________________________________________________________________________________
Mathematica [A]
time = 0.40, size = 149, normalized size = 1.77 \begin {gather*} \frac {\sqrt {2} \cosh (x) \sqrt {a+b \text {sech}^2(x)} \sqrt {a+b+a \sinh ^2(x)} \left (3 \sqrt {a} (a+b) \sinh ^{-1}\left (\frac {\sqrt {a} \sinh (x)}{\sqrt {a+b}}\right )-\sqrt {a+b} \text {csch}(x) \left (4 a+3 b+(a+b) \text {csch}^2(x)\right ) \sqrt {\frac {a+b+a \sinh ^2(x)}{a+b}}\right )}{3 (a+b)^{3/2} \sqrt {a+2 b+a \cosh (2 x)} \sqrt {\frac {a+b+a \sinh ^2(x)}{a+b}}} \end {gather*}
Antiderivative was successfully verified.
[In]
Integrate[Coth[x]^4*Sqrt[a + b*Sech[x]^2],x]
[Out]
(Sqrt[2]*Cosh[x]*Sqrt[a + b*Sech[x]^2]*Sqrt[a + b + a*Sinh[x]^2]*(3*Sqrt[a]*(a + b)*ArcSinh[(Sqrt[a]*Sinh[x])/
Sqrt[a + b]] - Sqrt[a + b]*Csch[x]*(4*a + 3*b + (a + b)*Csch[x]^2)*Sqrt[(a + b + a*Sinh[x]^2)/(a + b)]))/(3*(a
+ b)^(3/2)*Sqrt[a + 2*b + a*Cosh[2*x]]*Sqrt[(a + b + a*Sinh[x]^2)/(a + b)])
________________________________________________________________________________________
Maple [F]
time = 1.70, size = 0, normalized size = 0.00 \[\int \left (\coth ^{4}\left (x \right )\right ) \sqrt {a +b \mathrm {sech}\left (x \right )^{2}}\, dx\]
Verification of antiderivative is not currently implemented for this CAS.
[In]
int(coth(x)^4*(a+b*sech(x)^2)^(1/2),x)
[Out]
int(coth(x)^4*(a+b*sech(x)^2)^(1/2),x)
________________________________________________________________________________________
Maxima [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Failed to integrate} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(coth(x)^4*(a+b*sech(x)^2)^(1/2),x, algorithm="maxima")
[Out]
integrate(sqrt(b*sech(x)^2 + a)*coth(x)^4, x)
________________________________________________________________________________________
Fricas [B] Leaf count of result is larger than twice the leaf count of optimal. 924 vs.
\(2 (70) = 140\).
time = 0.49, size = 2341, normalized size = 27.87 \begin {gather*} \text {Too large to display} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(coth(x)^4*(a+b*sech(x)^2)^(1/2),x, algorithm="fricas")
[Out]
[1/12*(3*((a + b)*cosh(x)^6 + 6*(a + b)*cosh(x)*sinh(x)^5 + (a + b)*sinh(x)^6 - 3*(a + b)*cosh(x)^4 + 3*(5*(a
+ b)*cosh(x)^2 - a - b)*sinh(x)^4 + 4*(5*(a + b)*cosh(x)^3 - 3*(a + b)*cosh(x))*sinh(x)^3 + 3*(a + b)*cosh(x)^
2 + 3*(5*(a + b)*cosh(x)^4 - 6*(a + b)*cosh(x)^2 + a + b)*sinh(x)^2 + 6*((a + b)*cosh(x)^5 - 2*(a + b)*cosh(x)
^3 + (a + b)*cosh(x))*sinh(x) - a - b)*sqrt(a)*log((a*b^2*cosh(x)^8 + 8*a*b^2*cosh(x)*sinh(x)^7 + a*b^2*sinh(x
)^8 - 2*(a*b^2 - b^3)*cosh(x)^6 + 2*(14*a*b^2*cosh(x)^2 - a*b^2 + b^3)*sinh(x)^6 + 4*(14*a*b^2*cosh(x)^3 - 3*(
a*b^2 - b^3)*cosh(x))*sinh(x)^5 + (a^3 + 4*a^2*b + 9*a*b^2)*cosh(x)^4 + (70*a*b^2*cosh(x)^4 + a^3 + 4*a^2*b +
9*a*b^2 - 30*(a*b^2 - b^3)*cosh(x)^2)*sinh(x)^4 + 4*(14*a*b^2*cosh(x)^5 - 10*(a*b^2 - b^3)*cosh(x)^3 + (a^3 +
4*a^2*b + 9*a*b^2)*cosh(x))*sinh(x)^3 + a^3 + 2*(a^3 + 3*a^2*b)*cosh(x)^2 + 2*(14*a*b^2*cosh(x)^6 - 15*(a*b^2
- b^3)*cosh(x)^4 + a^3 + 3*a^2*b + 3*(a^3 + 4*a^2*b + 9*a*b^2)*cosh(x)^2)*sinh(x)^2 + sqrt(2)*(b^2*cosh(x)^6 +
6*b^2*cosh(x)*sinh(x)^5 + b^2*sinh(x)^6 - 3*b^2*cosh(x)^4 + 3*(5*b^2*cosh(x)^2 - b^2)*sinh(x)^4 + 4*(5*b^2*co
sh(x)^3 - 3*b^2*cosh(x))*sinh(x)^3 - (a^2 + 4*a*b)*cosh(x)^2 + (15*b^2*cosh(x)^4 - 18*b^2*cosh(x)^2 - a^2 - 4*
a*b)*sinh(x)^2 - a^2 + 2*(3*b^2*cosh(x)^5 - 6*b^2*cosh(x)^3 - (a^2 + 4*a*b)*cosh(x))*sinh(x))*sqrt(a)*sqrt((a*
cosh(x)^2 + a*sinh(x)^2 + a + 2*b)/(cosh(x)^2 - 2*cosh(x)*sinh(x) + sinh(x)^2)) + 4*(2*a*b^2*cosh(x)^7 - 3*(a*
b^2 - b^3)*cosh(x)^5 + (a^3 + 4*a^2*b + 9*a*b^2)*cosh(x)^3 + (a^3 + 3*a^2*b)*cosh(x))*sinh(x))/(cosh(x)^6 + 6*
cosh(x)^5*sinh(x) + 15*cosh(x)^4*sinh(x)^2 + 20*cosh(x)^3*sinh(x)^3 + 15*cosh(x)^2*sinh(x)^4 + 6*cosh(x)*sinh(
x)^5 + sinh(x)^6)) + 3*((a + b)*cosh(x)^6 + 6*(a + b)*cosh(x)*sinh(x)^5 + (a + b)*sinh(x)^6 - 3*(a + b)*cosh(x
)^4 + 3*(5*(a + b)*cosh(x)^2 - a - b)*sinh(x)^4 + 4*(5*(a + b)*cosh(x)^3 - 3*(a + b)*cosh(x))*sinh(x)^3 + 3*(a
+ b)*cosh(x)^2 + 3*(5*(a + b)*cosh(x)^4 - 6*(a + b)*cosh(x)^2 + a + b)*sinh(x)^2 + 6*((a + b)*cosh(x)^5 - 2*(
a + b)*cosh(x)^3 + (a + b)*cosh(x))*sinh(x) - a - b)*sqrt(a)*log(-(a*cosh(x)^4 + 4*a*cosh(x)*sinh(x)^3 + a*sin
h(x)^4 + 2*(a + b)*cosh(x)^2 + 2*(3*a*cosh(x)^2 + a + b)*sinh(x)^2 + sqrt(2)*(cosh(x)^2 + 2*cosh(x)*sinh(x) +
sinh(x)^2 + 1)*sqrt(a)*sqrt((a*cosh(x)^2 + a*sinh(x)^2 + a + 2*b)/(cosh(x)^2 - 2*cosh(x)*sinh(x) + sinh(x)^2))
+ 4*(a*cosh(x)^3 + (a + b)*cosh(x))*sinh(x) + a)/(cosh(x)^2 + 2*cosh(x)*sinh(x) + sinh(x)^2)) - 4*sqrt(2)*((4
*a + 3*b)*cosh(x)^4 + 4*(4*a + 3*b)*cosh(x)*sinh(x)^3 + (4*a + 3*b)*sinh(x)^4 - 2*(2*a + b)*cosh(x)^2 + 2*(3*(
4*a + 3*b)*cosh(x)^2 - 2*a - b)*sinh(x)^2 + 4*((4*a + 3*b)*cosh(x)^3 - (2*a + b)*cosh(x))*sinh(x) + 4*a + 3*b)
*sqrt((a*cosh(x)^2 + a*sinh(x)^2 + a + 2*b)/(cosh(x)^2 - 2*cosh(x)*sinh(x) + sinh(x)^2)))/((a + b)*cosh(x)^6 +
6*(a + b)*cosh(x)*sinh(x)^5 + (a + b)*sinh(x)^6 - 3*(a + b)*cosh(x)^4 + 3*(5*(a + b)*cosh(x)^2 - a - b)*sinh(
x)^4 + 4*(5*(a + b)*cosh(x)^3 - 3*(a + b)*cosh(x))*sinh(x)^3 + 3*(a + b)*cosh(x)^2 + 3*(5*(a + b)*cosh(x)^4 -
6*(a + b)*cosh(x)^2 + a + b)*sinh(x)^2 + 6*((a + b)*cosh(x)^5 - 2*(a + b)*cosh(x)^3 + (a + b)*cosh(x))*sinh(x)
- a - b), -1/6*(3*((a + b)*cosh(x)^6 + 6*(a + b)*cosh(x)*sinh(x)^5 + (a + b)*sinh(x)^6 - 3*(a + b)*cosh(x)^4
+ 3*(5*(a + b)*cosh(x)^2 - a - b)*sinh(x)^4 + 4*(5*(a + b)*cosh(x)^3 - 3*(a + b)*cosh(x))*sinh(x)^3 + 3*(a + b
)*cosh(x)^2 + 3*(5*(a + b)*cosh(x)^4 - 6*(a + b)*cosh(x)^2 + a + b)*sinh(x)^2 + 6*((a + b)*cosh(x)^5 - 2*(a +
b)*cosh(x)^3 + (a + b)*cosh(x))*sinh(x) - a - b)*sqrt(-a)*arctan(sqrt(2)*(b*cosh(x)^2 + 2*b*cosh(x)*sinh(x) +
b*sinh(x)^2 + a)*sqrt(-a)*sqrt((a*cosh(x)^2 + a*sinh(x)^2 + a + 2*b)/(cosh(x)^2 - 2*cosh(x)*sinh(x) + sinh(x)^
2))/(a*b*cosh(x)^4 + 4*a*b*cosh(x)*sinh(x)^3 + a*b*sinh(x)^4 - (a^2 + 3*a*b)*cosh(x)^2 + (6*a*b*cosh(x)^2 - a^
2 - 3*a*b)*sinh(x)^2 - a^2 + 2*(2*a*b*cosh(x)^3 - (a^2 + 3*a*b)*cosh(x))*sinh(x))) + 3*((a + b)*cosh(x)^6 + 6*
(a + b)*cosh(x)*sinh(x)^5 + (a + b)*sinh(x)^6 - 3*(a + b)*cosh(x)^4 + 3*(5*(a + b)*cosh(x)^2 - a - b)*sinh(x)^
4 + 4*(5*(a + b)*cosh(x)^3 - 3*(a + b)*cosh(x))*sinh(x)^3 + 3*(a + b)*cosh(x)^2 + 3*(5*(a + b)*cosh(x)^4 - 6*(
a + b)*cosh(x)^2 + a + b)*sinh(x)^2 + 6*((a + b)*cosh(x)^5 - 2*(a + b)*cosh(x)^3 + (a + b)*cosh(x))*sinh(x) -
a - b)*sqrt(-a)*arctan(sqrt(2)*(cosh(x)^2 + 2*cosh(x)*sinh(x) + sinh(x)^2 + 1)*sqrt(-a)*sqrt((a*cosh(x)^2 + a*
sinh(x)^2 + a + 2*b)/(cosh(x)^2 - 2*cosh(x)*sinh(x) + sinh(x)^2))/(a*cosh(x)^4 + 4*a*cosh(x)*sinh(x)^3 + a*sin
h(x)^4 + 2*(a + 2*b)*cosh(x)^2 + 2*(3*a*cosh(x)^2 + a + 2*b)*sinh(x)^2 + 4*(a*cosh(x)^3 + (a + 2*b)*cosh(x))*s
inh(x) + a)) + 2*sqrt(2)*((4*a + 3*b)*cosh(x)^4 + 4*(4*a + 3*b)*cosh(x)*sinh(x)^3 + (4*a + 3*b)*sinh(x)^4 - 2*
(2*a + b)*cosh(x)^2 + 2*(3*(4*a + 3*b)*cosh(x)^2 - 2*a - b)*sinh(x)^2 + 4*((4*a + 3*b)*cosh(x)^3 - (2*a + b)*c
osh(x))*sinh(x) + 4*a + 3*b)*sqrt((a*cosh(x)^2 + a*sinh(x)^2 + a + 2*b)/(cosh(x)^2 - 2*cosh(x)*sinh(x) + sinh(
x)^2)))/((a + b)*cosh(x)^6 + 6*(a + b)*cosh(x)*sinh(x)^5 + (a + b)*sinh(x)^6 - 3*(a + b)*cosh(x)^4 + 3*(5*(a +
b)*cosh(x)^2 - a - b)*sinh(x)^4 + 4*(5*(a + b)...
________________________________________________________________________________________
Sympy [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \int \sqrt {a + b \operatorname {sech}^{2}{\left (x \right )}} \coth ^{4}{\left (x \right )}\, dx \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(coth(x)**4*(a+b*sech(x)**2)**(1/2),x)
[Out]
Integral(sqrt(a + b*sech(x)**2)*coth(x)**4, x)
________________________________________________________________________________________
Giac [F(-2)]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Exception raised: TypeError} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(coth(x)^4*(a+b*sech(x)^2)^(1/2),x, algorithm="giac")
[Out]
Exception raised: TypeError >> An error occurred running a Giac command:INPUT:sage2:=int(sage0,sageVARx):;OUTP
UT:Error: Bad Argument Type
________________________________________________________________________________________
Mupad [F]
time = 0.00, size = -1, normalized size = -0.01 \begin {gather*} \int {\mathrm {coth}\left (x\right )}^4\,\sqrt {a+\frac {b}{{\mathrm {cosh}\left (x\right )}^2}} \,d x \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
int(coth(x)^4*(a + b/cosh(x)^2)^(1/2),x)
[Out]
int(coth(x)^4*(a + b/cosh(x)^2)^(1/2), x)
________________________________________________________________________________________